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^ ! Abstract 

In this paper, we calculate analytically the baryonic density and suscepti- 
bilities, which are sensitive probes to the fermionic degrees of freedom, in a 
^! holographic model of QCD both in its hot QGP phase and in its cold dense 

r-| \ phase. Interesting patterns due to strong coupling dynamics will be shown 

and valuable lessons for QCD will be discussed. 
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1. Introduction 



Historically string theory was developed during the attempt to build dual 
models for the strong interaction describing various hadrons and their scat- 
tering amplitudes [H- After the advent of Quantum Chromodynamics (QCD) 
as the fundamental gauge theory of strong interaction, string theory departed 
into a quite different route with many fascinating discoveries. More recently, 
a new twist has been added to the story between string theory and QCD: 
the AdS/CFT correspondence [2I, S, S], or more generally the gauge/gravity 
duality, has provided a way to study the strong coupling dynamics of non- 
Abelian gauge theories (like QCD) via its string theory dual in the much 
more tractable supergravity regime. A lot of interesting results have been 
achieved along this direction with applications in various aspects, for reviews 
see e.g. 0, 0, 0, i, S 0, [111, ^2 . 



While a specific gravity dual for real QCD hasn't been found yet, many 
important insights into the nonperturbative aspects of QCD have been ob- 
tained via gauge/gravity duality. Let's just mention two remarkable examples 
of such kind. Study of the so-called quark-gluon plasma (QGP), the decon- 
fined phase of QCD at temperature higher than the deconfinement transition 
T > Tc ~ 170 — 200MeV and low baryonic density, has been very impor- 
tant but proved to be rather difficult. Two powerful tools for such study 
include the lattice QCD and the heavy ion collisions experiments (e.g. at 
RHIC facility): interestingly, the former approach found that the the QGP 
thermodynamics (e.g. the energy density) deviates from Stefan-Boltzmann 
limit constantly by about 20% in the range 1.5 — 4Tc [l^, while the latter ap- 
proach found that the QGP produced at RHIC (corresponding to 1 — 2Tc) has 
an extremely small shear viscosity to entropy density ratio rj/s of the order 



0.1 [Ij]. Both discoveries challenged naive expectation of a weakly coupled 
QGP right above Tr, and led to the paradigm shift to a strongly coupled QGP 
(sQGP) 0, [lil . [Taj . It turns out that both results can be much better un- 
derstood in light of the strong couphng results from gauge/gravity duality: 
the thermodynamics of CFT plasma at infinitely strong coupling has been 
calculated via its dual black hole to be exactly | of the Stefan-Boltzmann 



limit [l7j, while a large class of strong coupling gauge theories with gravity 



dual have been shown to have their 2 = -L ~ g.OS (iSi. These successes have 

in 



inspired a lot of activity to extract useful information for QCD and QGP via 
gauge/gravity duality, see e.g. reviews in 0,[ii|,[i5j. 



In this paper, we focus on the fermionic sector i.e. the degrees of freedom 
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carrying baryonic charges in QCD thermodynamics, and use the gauge/gravity 
duahty to quahtatively obtain some of their nonperturbative features. There 
are very rich dynamics associated with the quarks in QCD, for example the 
spontaneous chiral symmetry breaking in QCD vacuum and its restoration 
at high temperature/density. There are also interesting phase structures in 
the low T and high baryonic density region of QCD phase diagram where the 
fermionic sector becomes dominant. Color superconductivity 19| was found 



to occur at very high density with the phase presumably a color-flavor-locking 



one 



201 ]. At moderately high density there are interesting phenomena like 



e.g. the interplay between the color superconductivity and chiral restora- 
tion and the pairing with mismatched Fermi surfaces 2l| 22] . More recently 
based on wisdom from large Nc argument, it has been proposed that there 
could be a new phase in the cold dense region, named "quarkyonic" phase 
231 ] 241 ]. which is confined and yet has thermodynamics scaling as Nc like a 
system made of quarks. There are also many recent results on the baryonic 
density and susceptibilities from lattice QCD which show nonperturbative 
patterns H S, H, [ij. 

It is thus of great interest to study these aspects of QCD with the handy 
tool of gauge /gravity duality. To this end, we use the Sakai-Sugimoto(SS) 



model 3l| , which have reproduced (at lease qualitatively) an impressive num- 
ber of QCD results, for example: the baryon property, form factor 3^], and 



the nuclear force 
ature 



331 ] ■ etc. The SS model has been studied at finite temper 



341 ] and at finite baryon density in [35|, |36|, |37j. We will particularly 



calculate the baryonic density and susceptibilities, which are sensitive probes 
to the fermionic degrees of freedom, from the Sakai-Sugimoto model both in 
the hot QGP phase and in the cold dense phase with the motivation to un- 
derstand their patterns under strong coupling and find valuable lessons for 
QCD. 

The paper is organized as follows. We give in Section. II a brief review of 
baryonic density and susceptibilities in QCD and in Section. Ill a brief review 
of pertinent Sakai-Sugimoto results obtained before. The baryonic density 
and susceptibilities in Sakai-Sugimoto model will be analytically calculated 
for the hot QGP phase in Section. IV and for the cold dense phase in Sec- 
tion. V. Finally in Section. VI we summarize the results and discuss relevant 
lessons for QCD. 
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2. Brief Review of Bciryonic Density and Susceptibilities in QCD 



2.1. Definition and Examples 

We start with a Taylor expansion of pressure P{T,n) with respect to 
chemical potential /i (with the convention that quark carries unit baryonic 
charge) at fixed T: 

P(T,,)^T±'-fi{!iy (1) 

n=0 

with the (dimensionless) baryonic susceptibilities dn{T) defined as 

Note for the above the odd-n ones vanish by symmetry. Furthermore 
we see all non-zero dn except n = represent certain contribution from 
the baryonic degrees of freedom, and importantly all non-baryonic degrees 
of freedom (e.g. the ghionic sector in QCD) do not directly contribute to 
them. So these derivatives probe the properties of the effective fermions in 
the system directly and sensitively. This can be seen also from the e.g. the 
baryonic density as given by 

To give an idea and to provide a benchmark of the baryonic density and 
susceptibilities, we explicitly evaluate these for a free gas of particles with 
mass M and baryonic charge B. The density is given by 



free 



rp3 PCX) C 



x^ x'^ 



(4) 

with P — 1/T and Ni denoting the number of internal degrees of freedom 
(e.g. color, flavor, spin, etc). The susceptibilities dn can be obtained by sub- 
sequent differentiation with respect to /i. We are particularly interested in 
two limiting cases: the non-relativistic(NR) limit with (3M oo and the 
ultra-relativistic(UR) limit with /3M — > 0. In these limits we can obtain 



3 



concrete results below. 

(i) NR limit: the baryonic density and susceptibilities in this limit are: 



free 



NR 



MT\^ _M 
27r ) 



e T 



Bfj. 
T 



NR 



4—) 

\27iT 



M' 



(5) 



(6) 



We observe a few important points at this limit: (a) all susceptibilities dn are 
positive and have the same dependence on T up to a constant coefficient ; (b) 
the ratio between successive susceptibilities is directly related to the baryonic 
charge carried by the degree of freedom, i.e. dn+2/dn = B^, independent of 
M,T and n; (c) for multi-component non-interacting gas of species Mi,Bi, 
all dn are simply a sum over species with the same formulae above and they 
remain all positive, but one then expect the ratios dn+2/dn —< Bf > to be an 
abundance-averaged baryonic charge which now depends on Bi and Mi,T,n 
as well. 

(ii) UR limit: the baryonic density and susceptibilities in this limit are: 



n 



UR 



free 



UR 



67r2 



B' 



+ 7r'B' 



?1 



di''^' 



B' 



UR 



jfree 
"'n>4 







UR 



(7) 



(8) 



Again for multi-component non- interacting gas of species M^, B^, the dn are 
simply a sum over species with the same formulae above. For such a Stefan- 
Boltzmann gas of quarks (with Bg — 1 convention) with spin Ng, flavor Nf 
and color N^, the susceptibilities per degrees of freedom (D.o.F) are simply 
dl^/{N,NfNc) = |, df^/{NsNfNc) = with all higher ones vanishing. 

In both limits of the free gas example above, we find all non-vanishing 
susceptibilities to be positive and proportional to which implies the con- 
tribution of fermions with large B becomes larger and larger with increasing 
order n. 

We end this part by emphasizing again that the density and susceptibil- 
ities are direct and sensitive probes to the fermions with baryonic charges 
in the system. In particular their deviation from the free patterns encodes 
important information about the dynamics and it is certainly of great inter- 
est to know the behavior of baryonic density and susceptibilities under strong 
coupling. 
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Figure 1: (from left to right) the susceptibilities ^2, d4, and de per D.o.F. calculated from 
lattice QCD. The dashed red lines in all panels are earlier 2-flavor results in [26|] while the 



solid blue lines in d2 and di plots are recent 3-flavor results by the same group in [27 1. 
The dotted green bands represent the Stefan-Boltzmann limit. 



2.2. The Susceptibilities in Lattice QCD 

While the susceptibihties contain very useful information about the fermionic 
degrees of freedom, it is generally hard to be calculated in a strongly coupled 
theory like QCD. Nevertheless such susceptibihties of QCD can be studied 
in its lattice formulation. In particular, calculating the dn (defined at zero 
/i) offers an important method (via Taylor expansion) to explore the QCD 
phase diagram at finite density which was traditionally unaccessible for lat- 
tice QCD due to the well-known "sign problem". These susceptibilities are 
also directly related to baryonic charge fluctuations in the thermal QCD mat- 



ter created in heavy ion collisions which can be experimentally measured [25 



The first lattice results for 2-flavor QCD 26|| with a relatively large pion mass 
by Karsch et al from a few years ago showed highly nontrivial patterns near 
Tc, with strong deviation from the free case in the region T ~ 1 — 2Tc. More 
recent 3-fiavor results with a more realistic pion mass by the same group 27 
still preserved the nontrivial patterns close to Tc though with deviation from 
the free case limited to be below ~ l.ATc. See also other lattice works in e.g. 



28|[29|[30| which all show similar behavior of the susceptibilities. 

We now make a detailed discussion on these lattice results. Let's start 
with the T < Tc part, i.e. the confined hadronic phase: in this phase we 
actually see from Fig{T] the same signs and similar T-dependence of ^2,4,6 
(including exponential-like growth close to T^) hinting at the NR limit of a 
free gas as in Eq.([6]). And indeed a simple hadronic resonance gas model in- 
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eluding various baryons from the Partiele Data Book ean nieely fit the lattice 



data for T < T^, see e.g. Fig.l in 



In contrast, the data for T > do 



not resemble any free gas (NR or UR) model at all. A few nontrivial features 
are readily visible from FiglH especially for 1 — lAT^. (i2,4,6 have rather dis- 
tinctive patterns, with d2 positive and mildly growing, d^ positive and rapidly 
dropping, while d^ negative and dropping even more abruptly. The suscep- 
tibilities results clearly disfavor any weakly interacting quasiparticle models 
for the deconfined quark-gluon plasma (QGP) phase right above Tc, instead 
they indicate non-perturbative effect from strong coupling. This observation 
turns out to be consistent with the conclusion from other independent ap- 
proaches that the QGP in the same 1 — 2Tc-region is rather strongly coupled, 
now called sQGP. In particular the experimental study of the QGP phase by 
heavy ion collisions at RHIC suggests that it behaves as a very good liquid 
Lattice works on e.g. thermodynamics 39|] and transport 
heavy quark potentials 41 1 and charmonium above Tc [4/ 



14, 15, 16 
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properties 

also point to a similar conclusion. A microscopic explanation of such strong 
coupling results may rely on understanding of the specific mechanism of how 
the confinement occurs toward from above: it has been suggested under 
the generic spirit of electric-magnetic duality that the QCD plasma close to 
Tc is actually a magnetic one, made of light and abundant monopoles which 



43, 44, 45, 46 



1.4Tc, a theoretical un- 



become Bose-condensed at and enforce confinement 

Returning to the non-trivial susceptibilities in 1 
derstanding is still lacking. One particular suggestion involves possible bound 
states of quarks and gluons even in the plasma phase due to the presence 



of still strong coupling in 1 — 2Tc 47|, |48|, |49|] , including both conventional 



colorless states (remnants of mesons and baryons) and colored states like di- 
quarks, q-g states and even polymer-like chains. Those bound states carrying 
baryonic charges can contribute substantially to the susceptibilities, as first 
pointed out by Liao and Shuryak in and especially become more and 
more important in higher orders. To see this, one just notes that diquarks 
have baryonic charge B = 2 and baryons B = 3, and that the susceptibilities 
go like dn ~ -B" as evident from Eq. (l^|5]) : this means even the density of the 
bound states may be small compared to quarks but they still can dominate 
e.g. d^, de and even higher order ones. It was shown in [sl] that the bound 
states contribution, mainly from baryons, could be dominant for the peak 
and wiggle structures seen in the lattice data. This was confirmed by an 
even better agreement with data in later studies using the so-called PNJL 
model 5^] in which due to Polyakov line suppression of colored states below 
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and around T^, precisely the 3-quark colorless combinations (e.g. baryons) 
dominate the nontrivial susceptibilities around Tg. Despite these progresses, 
it is still unclear how the strong coupling dynamics affects the susceptibilities 
from individual quarks — this is what we attempt to address in the present 
paper. 



To conclude this section, we have seen that the baryonic density and 
susceptibilities are sensitive probes to the degrees of freedom carrying bary- 
onic charges, and can exhibit rather nontrivial structures in strong coupling 
regime as shown by lattice QCD data. It is natural, then, to study these 
properties for strongly coupled gauge theories in general and for useful mod- 



els of QCD in particular. To this end, the Sakai-Sugimoto model [31 



IS a 



good choice: as a holographic model of QCD, it is calculable in strong cou- 
pling regime by virtue of the gauge/string duality and has been shown to 
have many realistic features of QCD. In the rest of the paper, we will cal- 
culate the baryonic density and susceptibilities in the Sakai-Sugimoto model 
and discuss the lessons for QCD. 



3. Brief Review of Sakai-Sugimoto Model with Chemical Potential 

In this section we summarize the Sakai-Sugimoto (SS) model for nota- 
tion and completeness. For a thorough presentation we refer 3l| for zero 



temperature and [3J] for finite temperature. 

The SS model, in brief, is defined by the dynamics of Nj D8-D8 branes 
in the background field (the metric, the dilaton, and the Ramond-Ramond 
field) generated by A^^^ D4-branes. In order not to disturb the background 
we require Nf -C iVc, which is called the "probe" limit and corresponds to 
the quenched approximation. The low energy dynamics of D8-D8 branes are 
governed by the Dirac-Born-Infeld (DBI) action and the Chern-Simons (CS) 
action: 

^DBi = -Tg j d\ e-"^ tr ^/-det{gMN + 27ca'FMN) , (9) 

Scs = C.trF^ , (10) 

where (f) and C3 are the dilaton and the Ramond-Ramond field. The metric 
generated by D4 brane is encoded in the induced metric gMN on the D8-D8 
branes. Fmn = Qm^n — On Am — i [Am, A]^] (M, N = 0, 1, ■ ■ ■ , 8) are the 
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field strength tensor of the U{Nf) gauge fields on the D8-branes. Note the 
tr is taken over the A^j- dimensional flavor space. Tg = (^2tt}»i1 tension 
of the D8-brane. In the following subsections we will specify the ingredient 
fields of the DBI and CS action: the induced metric, dilaton, RR field, and 
the gauge field. 

3. 1 . Ingredient fields 

The induced metric on the D8 branes from the D4 branes background 



metric can be written as 31, 34 



dsls = 9oo{dX§f + gUdXf + QuudU^ + Qssd^l , (11) 

wher^ 

a ^ 1 , 7 ^ + i^^f^ fiU) , (Confined phase) 

1 fdXi\'^ fU^ ^ 



tt ^ f{U) , 7 ^ + y'olfj y^j ■ (Deconfined phase) 

The embedding information is encoded only in 7 and thereby guu- For 
definition of the warping factor f{U) in both confined and deconfined phases, 
see the appendix. 

The dilaton and RR field are given as 

^' = 9s{^"\ F,^dC, = ^-^e,, (12) 

where gs is the string coupling constant, ^4 = 87r^/3 is the volume of the 
unit 5*^ and €4 is the corresponding volume form. 

For the gauge field we only consider the time component of the U{1) part 
of the U{Nf) gauge field, which is normalized as 

^MU) . (13) 



In the appendix, the background metric and coordinate notation convention are shown. 
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where we also assumed the Aq is only a function of U. This choice of the bulk 
gauge field is a standard holographic way to introduce the chemical potential 
in the boundary field theory. y4o(oo) of the classical solution is identified with 
the chemical potential. 

3.2. The DBI action 

the DBI action ([9]) reads 



DBI 



dUe 



^9ss9xx 



9oo9uu 



27ia'A, 



(14) 



where Ns = 2 and reflects the two contributions from D8 and D8 branes and 
^4 results from the trivial angle integration of S^. The symbol ' denotes the 
differentiation with respect to U. Note that googuu is the only place where 
the confined and deconfined phases are distinguished. 



In terms of the dimensionless variables defined as 36 



U 



the DBI action reads 



X4 



X4 

R 



R ' 



27ra'A 



'2NfR 



(15) 



5i 



DBI 



M / duu\ f\u){x'^{u)Y + 



1 



f{U) 

1 



confined 
deconfined 



where 



^^^^ _ NsNfTsR^Q^{(3V^) 
9s 



(16) 



j3 is the inverse temperature, V3 is the volume of X space. 

For simplicity in this paper we only consider the antipodal configura- 
tion, in which the D8 branes and D8 branes are maximally separated at the 
boundary. This implies that, in the confined phase, the D8-D8 branes are 
connected ai U = Ukk, and, in the deconfined phase, D8 and D8 branes 
are parallel to each other extending to the black hole horizon. Since it also 
imphes that X4 is constant {x'^ = 0) the DBI action is reduced to 



S 



DBI 



M duu-^.n /(f) -(w) 



confined 
deconfined 



(17) 
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The antipodal configuration simplifies the problem since we only need to 
solve the one variable (ao) equation instead of the coupled equations of two 
variables (ao, X4). It also simplifies the phase diagram of holographic QCD 
since there is no "cusp" configuration studied in [36|]. Thus the deconfined 
phase at T > Tc and any fi corresponds to the hot QGP phase, |fl and the 
confined phase at T < and sufficiently large fi > fic corresponds to a cold 
dense phase made of baryonic matter while at small fi < fic the phase is the 
trivial vacuum configuration with zero baryonic density. The onset chemical 
potential, fic, is defined in (l26l) . The antipodal configuration will allow us to 
obtain analytic results for both the hot QGP phase and the cold dense phase, 
which shall remain qualitatively the same for large non-maximal separations 
of D8-D8 branes. 

3.3. Grand potential 

The grand potential is identified with the on-shell DBI action. For con- 
venience we will compute the rescaled grand potential defined as [f] 

fl{t,fi) = ■^5'dBI I on-shell , (18) 

where t is the dimensionless temperature and /x is the dimensionless chemical 
potential defined as 

t = T-R , fi = ao(oo) . (19) 

3.3.1. QGP phase 

Let us first consider the deconfined phase, 

^DBi = ^f J duui ^i-{d'^{u)y . (20) 

Since Sq is a cyclic coordinate, its conjugate momentum is conserved, which 
is defined up to a normalization as 

d = -—p , . = , = , (21) 

M5a'^{u) Vl - {a',{u)f 



^In general there could also be other phases corresponding to connected configurations 
in deconfined phase as shown in [36] . 

•^The grand potential in SS model has been computed in several papers 
with various notational conventions. In this paper we follow the convention in 
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and do is easily solved. Especially 

^ = ao(oo) = [ du f , (22) 

Jut V + " 

where a{uT) = for the regularity at the horizon. With this solution the 
grand potential reads 

^(i,/^)=/ du , (23) 

Jut V ^ + "(Ai) 

where d{^) is the function of ^ via fl22]) . 

5". 5".^. CoW (iense phase 

In the confined phase we can do the same analysis with 



^DBi = M / duui ^f{u) - (a'o(M))2 . (24) 



However in this case we need to consider the explicit source term. Contrary 
to the deconfined phase, D8 and D8 branes are connected at m = ukk- For a 
nontrivial So there must be a singularity at u = ukk and to take this account 
we consider the D4 branes wrapping S*^ as the baryon source. The source 
of a uniform distribution of D4 branes (~ d) can be explicitly introduced by 



the CS action. Referring to [36| for more detail, we simply quote the final 



result for the grand potential and chemical potential 



1 ^5 



^^^^ ' '^'^ ^/m v^^^^^+diw ' ^^^^ 



Ukk 



/i = ao(oo) = / du^==^=== + /ie • (26) 

J UKK ^/f{u)Vu^ + d^ 

where /ic = \ukk is the onset chemical potential of the cold dense phase. 

4. Baryonic Density and Susceptibilities in the QGP Phase 

At T > Tc the Sakai-Sugimoto model has a deconfined, chirally symmet- 
ric QGP phase for any /i. We make use of the analytic results for QGP phase 
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to analyze its baryonic degrees of freedom. Using Eq. fl23f22l) and subtract- 
ing out the vacuum part in the grand potential, we explicitly write out the 
following equations as our starting point at given T, /i: 



Pqgp [T,d{T,fi)] 



2t. 2 , . 5,1 2 ^ ^ /I 1 6 

-r^rfs + -UT{d' + 4)5 - -urd 2F1 ( -, -; -; 



(27) 



r.d.-...F,(-,-;-;-^ -,^0, r.^^i^ (28) 



In principle we need to solve d{T, n) from the constraint and then obtain the 
full P[T,fi] to calculate density and susceptibilities. However, we can make 
use of the chain rule to do the calculation without solving the constraint. In 
other words, we (temporarily) consider fi as being determined by T, d: 

^ = r.d.-„..F.(-,-;-;--^j (29) 

And we first evaluate the derivatives over d: 

dP _2uT(d'y /3 3 13 d'\ 

-dd-^y^j ^^^^10' 2' 10' "4 J 

dfi^2uT(f_y p. /^A 3 13 _o?^\ 

dd 3d Ut/ VlO'2'10' u^^J ^ ^ 

The density is then given by = = d, as it should. To determine its 

ultimate dependence on T, we will have to solve the constraint. 

We now calculate the susceptibilities dn as defined in Eq.([2]). To do 
that we introduce Pn{T,d) = with dn = t"""^ ■ Pn{d — ^ 0) (noting that 
— i> is equivalent to ;U — > 0). P„ can be evaluated order by order, using 
1 ~ ^ail/dd ■ "^^^ ^"cicil results to the order n — 10 are given below: 



■ A\ -2n+7 

J _ -n+7/-2 ,n-4 _ / 1 + 

"n=2,4,6,8,10 — Xn "^T ~ I ~Q~ ) 



n+7/2 _ / _ ) +-n+3 

3 3^ 3^ ■ 5 ■ 31 



2 ' 23.13 ' ^" 25 -132 -23 ' 

_ 3*^-52 •7-3011 _ 3^ ■ 52 ■ 7 ■ 24546787 

~ 27. 11- 133 -23 ' ~ ~ 29 ■ 11 ■ 134 ■ 232 ■ 43 ^ ' 
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These results show very distinctive patterns: (a) first of all we see simple 
power dependence on temperature but different orders depend on T very 
differently; (b) particularly, only ^2 has positive power of T-dependence, 
i.e. growing with T, while all higher order susceptibilities have negative 
power thus vanish in the T — > cxd limit; (c) furthermore, we notice there is 
an alternating sign pattern for n > 2, i.e. 6/43 are positive while (Iq^iq are 
negative; (d) finally we notice the ratios between successive susceptibilities 
are dn+2/dn ~ t~'^- As is evident from the above results, even in the QGP 
phase with quarks as the basic baryonic charge carriers, the strong interaction 
modifies the behavior significantly, resulting in non-perturbative patterns of 
these susceptibilities. 

At this point, it would be interesting to see the actual dependence of the 
susceptibilities on physical parameters by recovering all the dimensions (and 
neglected constants) of involved quantities, i.e. P, n, T. This can be done via 
the following: 



P ^ Px 



28 • 3 ■ vrs . Ml ^ ■ 



A 



(33) 



We also re-scale temperature by t = TR ^ 
Eventually we arrive at the results below: 



dr. 



X 



4 ■ vr ■ Mkk ■ 

L 2n -J L 27V 



Cny<NsXN,xNfX ( ^ 



n— 3 



(34) 



with = Xn- 3^""*^/ (tt ■ 2"^'^) (note the sign patterns of ^„ follow from Xn, i-e. 
1^2,4,8,. .. > while ^6,10,... < 0). First of all we notice the degrees of freedom 
counting reflects quark-like dependence, as baryonic or any other non-single 
combinations of quarks would result in different dependence on both Nc and 
Nf. On the other hand we see nontrivial power dependence on the coupling: 
again ^2 ~ AT grows both with A and T but all higher susceptibilities will 
be suppressed at very strong coupling. 

Finally let's discuss the asymptotic behavior of baryonic density at very 
large and very small density /chemical potential, i.e. fi 00 and fi 0. In 
these limits we can solve the constraint equation to obtain the density. 



While we only worked out the susceptibilities explicitly to the 10th order, it is plausible 
that the dependence on A and T is valid also for even higher orders. 
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(i) Dense Limit in QGP Phase: In the dense hmit /i — oo and thus d — >■ 
oo, we solve the constraint equation Eg. (1251) to leading order and obtain the 
pressure and baryonic density to be (after recovering physical dimensions): 

7 5 
pdense / ii \ 2 pdense 7 , / 11 \ 2 

'(35) 

We can also calculate the energy per particle 

- an-P 5 20 - 

E^- = -n = —EsB 36 

with Es.B. = <ifi/4: representing the energy per particle for a free gas of 
massless fermions at zero temperature and finite /i. Amusingly, the strong 
coupling result misses the non-interacting limit by only 1/21 ~ 5% despite 
the strong coupling dynamics. 

(ii) Dilute Limit in QGP Phase: In the dilute limit /i — and thus 

(i — s> 0, we note that 

d^O 

can then calculate the baryonic density to leading order 



Fi (ii;f;-f)][^^-r^-f -i4 + 6(rf^)and 



^dilute / T \ / II 

r^N.N.NAU,] (^) (37) 



This result has linear dependence on fi, in common with the leading order 
at small fi of free fermion gas at UR limit. The strong coupling density here 
however depends on coupling A and also has dependence, differing from 
the free case which has no coupling and has dependence. 



5. Baryonic Density and Susceptibilities in the Cold Dense Phase 

We now turn to the cold dense phase at T = 0. To study this phase, we 
focus on the situation with maximal D8-D8 separation which allows tractable 
analytic formulae and much simplifies the calculation. For non-maximal 
separation, the physics at T = remains qualitatively the same while the 
calculation is much more involved. 

We start with Eg. fl25ll26p and further rewrite the pressure and chemi- 
cal potential as a function of d into the following forms by re-scaling the 
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quantities with ukk'- 



[ukk] 



du 



\fuF 



1 + 4 



(38) 



u= = /ic + / g-u , = -\ / ^ (39) 

In the above (i^ = jv^xK ^^"^ A'-c = IJ'cIukk = 1/3- Also note in the pressure 
we have subtracted out the vacuum part, e.g. the (—1) within the bracket. 
These can be easily solved numerically: for each given jl we can first solve d 
from Eq. fl39p and then obtain P from Eq. (l38l) . and the results are shown as 
solid blue lines in Fig. 2. 

A few comments are in order. First we notice that for fi < fie there is 
no solution for d as is evident from Eq. fl39|) . so the solution is the trivial 
U-shape and the system is in vacuum phase without any dependence on jl. 
When ft > fie baryons start to emerge and a new phase with nonzero baryonic 
density takes over the vacuum one. Furthermore both the pressure and the 
baryonic density are continuous at fic transition while the first derivative of 
d versus /i (i.e. the lowest susceptibility at fic) is not, which implies a second 
order phase transition. In physical unit, one has /Xc = = ^ x A x T^: 

with A ~ 20 — 50, the relation reasonably agrees with current rough estimate 
of fic/Tc in QCD. 

We now make an expansion of the pressure and density in the cold dense 
phase close to fic- This can be done by systematically analyze the expansion 
of 0? — ^ order by order in the integrands of Eqs. (l38ll39p . The result for the 
pressure is: 

(40) 

r(ir ' r(f) _ 

We skip to show the higher order coefficients whose expressions become 
really long. Here we give an idea of the numbers: C2 = 0.477465, C4 = 



P = 2jc2n(/^-/^c)^" 

n 

_ 3 _ 9r(f) _ 243 
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Figure 2: Comparison of the full numerical results(solid blue line) for the pressure (left) and 
density(right) with their respective Taylor series results (dashed lines with growing lengths) 
truncated at the Ist(red), 2nd(green), 3rd(orange), 4th(magenta) , and 5th(purple) order. 
The vertical dashed black line indicates the position of the phase transition. 



0.0441772, cg = 0.001576, Cg = -0.000365984, cio = 0.0000107106. These co- 
efficients can be directly related to baryonic susceptibilities defined at /i = /Ic, 



I.e. 



D2n = 



92"(p//i: 



9(/i//Xe)2- 



~^2n-4 Q2np 



A'=A»c 



U 



1 
2 

KK 



u 



1_ 

2 

KK 



{2n)\c2r. 



(41) 



In Fig.2(left) we compare the full numerical result (solid blue line) of the 
pressure with its Taylor series results truncated at the Ist(red), 2nd(green), 
3rd(orange), 4th(magenta), and 5th(purple) order respectively. Interestingly 
we found that the expansion till the 2nd order, i.e. the series with {jl — /I^)^ 
and (yU — /ic)^ terms only, agrees remarkably well with the full result up to 
very large /i, indicating some fine cancellations among all higher order terms. 
The expansion for the density is given by: 



d 



(42) 



3 f 

— , /3 — ^ 
TT 7r2 



36 r(f) 



729 



2r(f)^ 0Fr(f) 
r(|)' r(f) 
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1.0 




0.6 ^, , ;, , , I , , , , I , , , , I , , , , I ] 

1 2 3 4 5 

ft 

Figure 3: The energy per particle versus chemical potential. The horizontal dashed green 
line indicates the free value. The vertical dashed black line indicates the position of the 
phase transition. 

Again we skip to show the higher order coefficients and give an idea of the 
numbers: /i = 0.95493, /a = 0.176709, /s = 0.00945598, fj = -0.00292787, /g = 
0.000107106. From the expansion it is clear that there is a jump in the slope 
of d versus /x, i.e. at /i — >■ /Ic~ and ^ at /i — fic^. We also note that these 
coefficients are related to those in the pressure expansion by /2n-i = (2n) ■C2n 
order by order as required hj d = dj^P. In Fig.2(right) we compare the full 
numerical result of the density (solid blue line) with its Taylor series re- 
sults truncated at the 1st (red), 2nd(green), 3rd(orange), 4th(magenta), and 
5th(purple) order respectively. Not surprisingly we again found that the ex- 
pansion till the 2nd order, i.e. the series with (/i — /i^) and (/2 — fic)^ terms 
only, agrees remarkably well with the full result up to very large ji. 

Finally we study another quantity of interest, i.e. the energy per particle 
E which at T = is given by ^ = {jld — P)/d = g{fi)fl. The function g{jl) 
is plotted in Fig. 3 as solid blue line. The dashed green line indicates the 
free case g = 3/4. We see that close to the transition point the energy per 
particle deviates much from the free value and drops very sharply, while at 
much larger density it curves back and approaches the free value from below 
very slowly. We notice that the curve crosses the free line at about ft ^ 2ftc: 
it may imply that below this density the system has repulsive interaction 
and above it the system has attractive interaction. In either regions, the 
deviation of E from free case is rather modest. 
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We also present the results with physical scales recovered as in Eq. 



X AT, X AT, X iV. X (A'A4,) X X: |J - 1 

4 / \2n 



2n 



(43) 



and 



D2n = 



a2"(P/^; 



where we have used fic = 



= NsXNfXN,x ^^^^^ X i (44) 
/.=Me ^ 32"-4(27r) A ^ ^ 



3 277r • 



6. Summary and Discussions 

In summary, we have calculated the baryonic density and susceptibilities 
in a holographic model of QCD. The results both for the hot QGP phase 
in Eq. (!34|) and for the cold dense phase in Eq. (l44|) show interesting patterns 
due to strong coupling dynamics. In the following we discuss relevant lessons 
for QCD that can be learned from the results. 

We first discuss the results for hot QGP phase close to T^. To give an 
idea of the numbers, we re-write the results in Eq. flMl) as: 



0.012 ■ AT 



N.NfN, 



0.37 



dp. 



26 



A3T3 



(45) 



We first note that a few qualitative features agree well with the lattice QCD 
data shown in FigH] the sign pattern of ^2.4,6 is the same, d2 approximately 
exhibits linear growth with T close to Tc while d^^ shows quick decrease 
with T, and also d^ vanishes much more abruptly than ^4. If one takes a 
A ~ 10 - 200 our d2 values are very close to the lattice results, but the 
(i4 and rfg values are too small. The comparison seems to indicate that 
dynamics from very strong coupling of quarks only tends to suppress higher 
order susceptibilities and may not be adequate to account for the d^ g close to 



- 17 from [31|. 
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Tc seen in lattice QCd[§ : this is yet another indication that those higher order 
susceptibihties may be attributed to multi-quark correlations like baryons 
which manifest themselves more and more in higher order susceptibilities as 
we already pointed out in Section. II. An even more critical quantity to further 
distinguish the two types of contribution to susceptibilities may be the non- 
diagonal susceptibilities 25| after introducing isospin chemical potential for 
quark flavors. These can be implemented in the holographic model of QCD 



5l[ |. and will be studied elsewhere. These susceptibilities can also be suitably 
combined to give the baryonic charge fluctuations in the QGP, which have 
phenomenological importance and can be measured in RHIC experiments 
[2 5( 1 : our results imply that strong coupling of quarks may enhance event-by- 
event fluctuation of baryon number as seen in the covariance < 5B'^ >~ AT 
while suppress even higher order variances. 

Another interesting observation is that for both phases we have noticed 
that in the dense regime the energy per particle E is very close to the Stefan- 
Boltzmann limit, i.e. (the so-called Bertsch number) is close to 1 de- 
spite the strong coupling. We remind that the result of entropy in A^^^ D3 
branes system (roughly involving adjoint gauge sector but without funda- 
mental fermions) deviates from the Stefan-Boltzmann limit by 25%. This 
comparison indicates that the strong coupling dynamics may somehow mod- 
ify the thermodynamics of fundamental fermions much more mildly than of 
the gauge fields. Interestingly lattice QCD seems to show similar situation: 
while for fermionic thermodynamics approaches the Stefan-Boltzmann 
limit already around 1.4Tc, the pure gauge lattice results for gluonic thermo- 
dynamics show no tendency toward the same limit even at 4Tc. 

Finally we point out that the cold dense phase of the Sakai-Sugimoto 
model studied above has its many features similar to the quarkyonic matter 



in large N^. QCD [23| [2J]. Both are in the confined regime but have their ther- 
modynamics scaling as NcNf like a quark system. The transition from the 
vacuum phase to the dense phase for both happens at about the baryon mass 
threshold 0, after which the baryonic density starts growing. While calcula- 



^As a caveat, one should be aware that already around the Tc scale, the many (un- 
included) higher KK modes may start kicking in and contributing to the susceptibilities — 
these become important with even higher temperature. At present it is not clear whether 
including their contributions would improve the agreement with lattice QCD results. 

'^The transition order however is the 1st in QCD while the 2nd in the present holo- 
graphic model: the difference may be due to the absence of nucleon-nucleon interaction 
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tion is hard for the quarkyonic matter in QCD, one may study its holographic 
dual quantitatively and obtain useful results as ours in Section. V. It will be 
very interesting to further investigate via the holographic method the many 
interesting features of the quarkyonic matter, e.g. the chiral symmetry and 
excitations near or deep into the Fermi surface, etc, which were previously 
qualitatively studied. 
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A. The Background Metric Generated by D4 branes 

In this appendix we summarize the metric generated by Nc D4 branes 
presented in [52 . 

At zero temperature, the metric generated by Nc D4-branes are given by 
ds' = (^^y {-{dXo)' + {dXf + f{U){dX,f^ 



where 



fiU) = 1 - ^ • (46) 

— * 

X = 2,3, while U{> Ukk) and ^4 are the radial coordinate and four angle 
variables in the Xs^ej.s.g direction. R is given by 

= TrgMl , (47) 



in the latter approach where the baryons are only introduced as static sources. See also 
comments in 13611. 
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where gs and Is are the string couphng and length respectively. Since X4 
is compactified to a circle, to avoid a conical singularity at U = Ukk, the 
period of 6X4^ is set to 

= • (48) 

Ukk 

The field theory is defined by Kaluza-Klein mass (Mkk) and the four-dimensional 
coupling constant at the compactification scale, gvM, 

O-rr 3 77^/2 2 

^ = 2^ ' = ^ = 27rMKK^?./s . (49) 

where 5'5(= {2'KYgsls) is the five-dimensional coupling constant obtained from 
D4 brane DBI action. From P7|) and (H9|) the parameters i?, J/kk; and gs 
may be expressed in terms of Mkk, ^(= 9ymXc), and as 



At finite-temperature, there are two possibilities [53|. One is to follow the 
standard prescription of the finite temperature field theory. The geometry is 
the same as the zero temperature apart from the fact that the time direction 
is Euclidean {Xq X^ = iXq) and compactified with a circumference = 
1/T: 

ds' = (^^y [(dX^)' + [dXf + frndX^f) 

+ T7 ^TTTTT + U^dnl 



u) \f{U) 

This corresponds to the confined phase which is thermodynamically preferred 
(i.e. with the smallest action) in the low temperature regime. The other 
possible geometry contains the black hole, which is another saddle point of 
the Euclidean path integral over supergravity configuration. The pertinent 
background is 

3/2 



ds^ = (/(t/)(rfXo^)^ + {dXy + {dX,f 

( du^ 
7 \fiu) 



+ 1^1 I ^TTTTT + U'dn^, 



21 




fiU) ^ 1 



(51) 



This corresponds to deconfined phase which is thermodynamically preferred 
(i.e. with the smallest action) at high temperature. To avoid a conical 
singularity at U = Ut the period of SXq of the compactified tE direction is 

set to 



which is identified with the inverse temperature. The confinement / deconfinement 
phase transition, i.e. the switching from one background geometry to the 
other, occurs when 6X4 = 6Xq i.e. at the critical temperature Tc, 
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